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SUMMARY 


I 

h 


At  millimeter  wave  frequencies  conventional  metal  waveguide  is  extremely 
lossy  and  very  expensive.  Dielectric  waveguide  may  be  an  inexpensive  and 
practical  substitute  for  many  applications. 

In  this  report,  we  investigate  the  properties  of  the  hybrid  HE^  mode  of  a 
cylindrical  dielectric  guide.  In  particular,  the  various  forms  of  loss  are  dis- 
cussed. Different  models  are  utilized  to  describe  loss  due  to  absorption  and 
scattering,  radiation  from  bends,  and  radiation  due  to  surface  imperfections. 

A heat  transfer  model  is  used  to  predict  the  power  handling  capabilities  of  a 
dielectric  rod. 

To  study  the  performance  of  an  actual  dielectric  guide,  an  experiment  was 
done  to  measure  the  dielectric  constant  and  attenuation  factor  of  a dielectric 
rod  propagating  a single  mode  at  70  GHz.  The  attenuation  factor  is  related  to 
an  equivalent  loss  tangent  that  takes  into  account  all  forms  of  loss.  By  com- 
paring this  loss  tangent  with  available  data  for  the  intrinsic  tan  5 of  the 
dielectric  material,  one  can  estimate  the  relative  importance  of  radiation 
losses  to  those  due  to  absorption, 

The  measurement  technique  utilizes  a long  section  of  dielectric  waveguide 
that  can  be  probed  anywhere  along  its  length.  The  standing  wavelength  along 
the  rod  is  easily  measured,  and  leads  to  a determination  of  the  dielectric 
constant.  Measurement  of  the  average  power  on  the  line  as  a function  of  position 
provides  the  attenuation  factor. 

Measurements  were  done  using  teflon,  polystyrene,  and  fused  quartz  rods. 
Dielectric  constant  at  71.0  GHz  was  measured  for  all  three  materials.  The 
attenuation  was  measured  for  several  different  polystyrene  and  teflon  rods. 

The  results  are  consistent  with  published  data  for  the  materials  at  microwave 
frequencies.  In  particular,  the  results  indicate  that  radiation  losses  may  be 
insignificant. 

The  report  concludes  with  a summary  of  the  theory  of  the  dielectric  wave- 
guide, and  a discussion  of  the  problems  that  remain  to  be  solved  before  dielectric 
waveguide  will  become  practical  for  use  with  millimeter  waves. 
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An  appendix  includes  numerical  calculations  of  the  propagation  vector  versus 
rod  diameter  for  different  values  of  e corresponding  to  teflon,  fused  quartz, 
and  polystyrene. 


Hi  Mm  y* 

PAUL  R.  WESSEL 
By  direction 
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Chapter  1 
INTRODUCTION 

The  use  of  millimeter  waves  for  radar,  radiometry,  and  communications  systems 
is  becoming  increasingly  popular.  Because  of  their  small  wavelength,  millimeter 
waves  provide  greater  resolution  for  tracking  purposes  than  do  lower  frequency 
microwaves.  There  is  no  shortage  of  possible  applications  of  millimeter  waves, 
but  there  are  a great  number  of  technical  problems  to  be  overcome  before  many  of 
these  applications  become  practical. 

Roughly  speaking,  millimeter  waves  fall  between  30  and  300  GHz.  Conventional 
microwave  techniques  can  be  used  at  the  lower  end  of  this  spectrum,  but  they  become 
impractical  as  the  frequency  increases.  In  particular,  conventional  metal  wave- 
guide becomes  extremely  lossy,  and  the  tolerances  required  for  manufacture  cause 
the  cost  to  be  prohibitive.  For  example,  the  theoretical  attenuation  of  a 
rectangular  guide  at  12  GHz  is  about  .03  dB/ft.  At  30  GHz,  it  is  about  .15  dB/ft, 
and  at  70  GHz  is  .5  dB/ft.  In  practice,  the  attenuation  is  about  twice  as  much 
due  to  mechanical  imperfections  and  contamination  of  the  conducting  surfaces  in 
the  guide.  There  is  one  notable  exception,  however.  The  TEg-j  mode  of  a 
cylindrical  guide  has  extremely  low  loss  due  to  the  absence  of  tangential  H 
field  at  the  walls  of  the  guide.  The  use  of  cylindrical  metal  guide  for  a long 
distance  communications  system  is  being  investigated  at  Bell  Laboratories.^  In 
general,  the  high  attenuation  makes  metal  waveguide  impractical  for  many  large 
scale  systems.  The  lack  of  high  power  sources  at  these  frequencies  only  compounds 
the  problem. 

There  are  several  alternatives  to  metal  waveguide  at  millimeter  wave  fre- 
quencies. They  include  microstrip,  strip-line,  image  line,  finline,  microguide. 


1. 


Abele,  T.  A.  et  al , "A  High-Capacity  Digital  Communication  System  Using  TEg, 
Transmission  in  Circular  Waveguide,"  IEEE  Trans.  Microwave  Theory  & Tech., 
Vol . MTT-23,  pp.  326-333,  April  1975. 
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and  dielectric  waveguide.  Dielectric  waveguide  isn't  nearly  as  sophisticated  as 
the  other  structures,  and  it  has  several  drawbacks.  Among  other  things,  it 
doesn't  interface  with  circuitry  easily,  it  can't  be  etched,  and  it  has  poor 
isolation  characteristics . 

However,  dielectric  waveguide  merits  study  because  it  is  so  simple  and 
inexpensive.  And,  in  theory,  its  attenuation  can  be  made  arbitrarily  small  by 
reducing  the  dimensions  of  the  guide. 

Dielectric  waveguide  has  been  the  subject  of  theoretical  analysis  since  1910, 

2 

when  it  was  investigated  by  Hondros  and  Debye.  In  recent  years,  it  has  been 
used  experimentally  for  measuring  dielectric  constants  of  materials  at  microwave 

3 

frequencies.  Dielectric  waveguide  will  soon  be  used  in  the  form  of  optical  fibers 
for  telephone  communications,  and  promises  to  revolutionize  the  communications 

4 

industry. 

In  this  technical  report,  we  investigate  the  properties  of  the  cylindrical 
dielectric  waveguide  at  millimeter  wave  frequencies.  The  objective  is  to  study 
the  problems  that  are  characteristic  of  any  waveguiding  structure  in  terms  of  the 
dielectric  cylinder,  so  that  one  may  develop  a basic  understanding  of  the 
advantages  and  limitations  of  any  type  of  structure.  The  dielectric  cylinder 
is  compared  with  metal  waveguide,  and  with  optical  fibers  when  relevant.  The 
discussion  is  limited  primarily  to  the  properties  of  the  HE^  mode  of  the  round 
dielectric  guide. 

The  technical  report  starts  with  a discussion  of  the  theories  of  the  lossless 
dielectric  slab  and  dielectric  cylinder.  The  general  solution  to  Maxwell's 
equations  is  derived  in  terms  of  a discrete  number  of  guided  modes,  and  an 
infinite  continuum  of  radiation  modes.  The  cutoff  and  propagation  characteristics 
of  the  guided  modes  are  discussed  in  terms  of  a guidance  condition.  Particular 
emphasis  is  placed  on  the  properties  of  the  dominant  HE^  mode  of  the  dielectric 
cylinder. 


2.  Hondros  and  Debye,  "Electromagnetische  Wellen  an  dielektrischen  Drahten," 
Ann.  der  Phys.,  Vol . 32,  pp.  465-476,  June  1910. 

3.  Cullen,  A.  L.,  and  Gillespie,  E.  F.  F.,  "A  New  Method  for  Dielectric 
Measurements  at  Millimeter  Wavelengths,"  Presented  at  the  Symposium  on 
Millimeter  Waves,  Polytechnic  Institute  of  Brooklyn,  April  1959. 

4.  TIME  Magazine,  Science  Section,  March  22,  1976. 
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The  effects  of  loss  are  introduced  using  a perturbation  technique.  The 

attenuation  in  a dielectric  rod  waveguide  due  to  absorption  and  scattering  is 

5 

described  using  the  theory  developed  by  Elsasser  *in  1949.  The  theory  of 
radiation  loss  from  bends  and  surface  imperfections  developed  by  Marcatili  and 
Marcuse  at  Bell  Laboratories  is  briefly  discussed. 

To  get  a practical  feel  for  the  problems  involved  in  using  a dielectric  guide, 
experimental  observations  of  a cylindrical  guide  were  made  at  70  GHz.  To  verify 
that  the  guide  was  performing  in  accordance  with  theory,  measurements  were  made 
of  the  dielectric  constant  and  attenuation  factor  of  several  teflon,  polystyrene, 
and  fused  quartz  rods.  By  relating  the  results  to  the  bulk  dielectric  constant 
and  loss  tangent  of  the  material  being  used,  one  can  verify  that  the  rod  was 
operating  in  the  HE^  mode.  In  addition,  the  experimental  determination  of  complex 
dielectric  constant  is  in  itself  an  important  result,  as  there  is  not  much  data 
available  for  materials  at  these  frequencies. 

Numerical  calculations  of  the  propagation  characteristics  of  the  HE^  mode 
are  included  in  an  appendix. 


! 


5.  Elsasser,  W.  M. , "Attenuation  in  a Dielectric  Circular  Rod,"  J.  Appl.  Phys.  20,  j 

p.  1193,  December  1949. 

6.  Marcatili,  E.  A.  J.,  "Bends  in  Optical  Dielectric  Guides,"  Bell  System  Technical 

Journal,  pp.  2103-2132,  September  1969.  I 

7.  Marcuse,  D. , and  Derosier,  R.M.,  "Mode  Conversion  Caused  by  Diameter  Changes 
of  a Round  Dielectric  Waveguide,"  B.S.T.J.  48,  No.  10,  pp.  3217-3232, 

December  1969. 

8.  Marcuse,  D. , "Radiation  Losses  of  the  Dominant  Mode  in  Round  Dielectric 
Waveguides,"  B.S.T.J.,  pp.  1665-1693,  October  1970. 

' 9.  Marcuse,  D. , "Radiation  Losses  of  Dielectric  Waveguides  in  Terms  of  the  Power 
Spectrum  of  the  Wall  Distortion  Function,"  B.S.T.J.,  pp.  3233-3242,  December 
1969. 
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Chapter  II 

THEORY  OF  THE  DIELECTRIC  WAVEGUIDE 

2.1  THE  DIELECTRIC  SLAB 

To  understand  the  behavior  of  a dielectric  waveguide,  it  is  necessary  to 
solve  Maxwell's  equations  subject  to  the  boundary  conditions  imposed  by  the 
presence  of  the  dielectric.  For  a non-conducting  dielectric  in  a region  where 
there  is  no  free  charge  or  current.  Maxwell's  equations  are 


\ 1 y x E + SR  = 0 

(2-la) 

V • B = 0 

(2-lc) 

at 

} j 

H 1 

v X H - — =0 

(2-lb) 

V • D = 0 

(2-1 d ) 

at 

where  D = eE 

B = uH 

(2-2) 

To  solve  the  field  equations  for  a dielectric  waveguide  of  uniform  cross- 
section,  it  is  convenient  to  use  the  approach  developed  by  Stratton.10  It  is 
possible  to  specify  both  the  electric  and  magnetic  fields  in  terms  of  a set  of 
Hertz  vectors  W and  n*,  such  that 


E = 

E 1 + E 2 

(2-3a) 

H = H 1 + H 2 

(2-3b) 

E1 

II 

< 

X 

< 

X 

(2-4a ) 

H1  = ev  x M 

(2-4c) 

at 

— * 

E2 

- -uv  X ^ 

(2-4b) 

-7  _ * 

N =vxvxn 

(2-4d) 

at 


10.  Stratton,  J.  A.,  Electromagnetic  Theory,  McGraw-Hill,  New  York,  1941. 


LL 
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_ _ * 

where  fl  and  n are  specified  by 

V2n  - ue  — =0  (2-5a)  V2n  * - ue  — =0  (2-5b) 

3t2  at-' 

The  vector  n can  be  related  to  the  conventional  scalar  and  vector  potentials 

by 

* = -V  • if  (2-6a)  A = -ue  — (2-6b ) 

at 

Equations  (2-3),  (2-4),  and  (2-5)  constitute  a general  solution  to  Maxwell's 
equations  in  homogeneous,  non-conducting  media  where  p and  J are  zero. 

To  describe  the  behavior  of  a dielectric  waveguide,  it  is  sufficient  to  find 
a set  of  Hertz  vectors  that  are  oriented  parallel  to  the  axis  of  the  guide.  Since 
these  vectors  are  unidirectional,  equations  (?-5)  become  scalar  equations.  As  an 
illustrative  example,  consider  the  slab  waveguide  of  Figure  2-1. 

We  wish  to  describe  a set  of  waves,  propagating  in  the  z direction,  that  are 

it 

guided  by  the  slab.  To  simplify  the  discussion,  we  will  treat  n and  n as  being 
dimensionless.  The  mathematics  will  be  unchanged,  and  appropriate  scaling  factors 
can  always  be  added  to  the  final  solution.  Let  n = n = c z.  Then  (2-5)  becomes 


Since  c is  different  in  the  two  regions  | x ) ^ d and  ) x | > d,  it  will  be 
necessary  to  solve  (2-7)  in  each  of  the  regions  using  the  proper  value  of  c. 


Since  we  are  interested  in  time-harmonic  waves  propagating  in  the  z direction, 
let 
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Expressing  equations  (2-4)  in  a Cartesian  system  and  making  use  of  (2-8), 
we  get 


<• 

±ih  !*■ 
dx 

= ±ih 

3y 

(k2 

- h2)<|< 

(2-9a) 

"1- 

ik2  dip 
yuj  3y 

ik2  dip 
yu)  3y 

"1  = 

0 

(2-9b) 

Ex- 

!y 

= -1va  u 

0 

(2-10a) 

H2 

y 

-*ihfj 

(k2 

- h2)4» 

(2-1  Ob) 

k = 

and 

ifi 

+ 

+ (k2 

- h2)*  * 

= 0 

(2-11) 

3X2 

3y2 

H 


and  E = E1  + E 2,  H = H1  + H 2 . 

The  significance  of  this  approach  can  be  seen  from  equations  (2-9),  (2-10), 
and  (2-11).  By  solving  a single  scalar  homogeneous  wave  equation,  we  are  able 
to  find  the  general  solution  for  the  electric  and  magnetic  fields  of  a dielectric 
waveguide.  Furthermore,  the  solutions  are  conveniently  separated  into  TE  and  TM 
waves.  (In  general,  there  is  no  TEM  solution  for  a dielectric  guide,  due  to  the 
impossibility  of  satisfying  all  of  the  boundary  conditions  simultaneously.) 

To  continue  the  analysis  of  the  dielectric  slab,  solve  equation  (2-11)  and 
impose  the  boundary  conditions  that  tangential  E and  H are  continuous  at  the 
surface  of  the  slab.  To  simplify  the  analysis,  let  us  restrict  the  solution  to 
be  transverse  electric  (TE)  and  to  be  independent  of  y (i.e.  - 3/3y  = 0).  Making 
use  of  equation  (2-8),  equation  (2-11)  becomes 


i 
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' . hJ)f(x)  = 0 (2-1 

9x2 


and  has  the  solutions 


s infix 

e"yx 

f ( X ) = 

COSBx 

eyx  _ 

(2-1 

where 

3 =^k2  - h2 

y =. 

^h2  - 

k2 

(2-1 

To  describe  a guided  mode. 

we  would 

like 

the  fields 

to  decay  with  increasing 

x,  so  that  there  is  no  electromagnetic  radiation  at  large  distances  from  the  slab 
Therefore  we  let 


Ae"' x 

x > d ■ 

f ( x ) = 

Ae*'* 

x < -d 

Bs  infix 

| x | < d_ 

fields  are 

x < -d 

1 x | < d 

x 

s d 

= -ipa)yAeyX 

II 

>> 

UJ 

-iiu.1  .BcosBx 

Ey 

iiiuiyAe  y X 

* - y 2 Ae  yx 

Hz  x 

B'BsinBx 

-y2Ae'yX 

= tih>Aeyx 

Hx  = 

♦ ihBBcosBx 

H = 

X 

Hh>  Ae*y  x 

where  exp  i'ihz  - i k>t > is  understood. 
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Matching  tangential  E and  H at  x = d, 


3Bcos8d  = -yAe  ' 
B*BsinSd  = -y’Ae  ^ 


(2-1 7a ) 
(2-1 7b) 


Dividing  ( 2- 1 7b ) by  (2-1 7a ) yields 

tamid  = y/8  (2-18) 

where  y and  8 are  functions  of  h defined  by  (2-14). 

Equation  (2-18)  is  the  guidance  condition  that  determines  the  eigenvalues  h, 
the  propagation  constant  in  the  z direction.  Since  the  tangent  function  is 
periodic,  there  are  numerous  values  of  h,  each  corresponding  to  a different  TE 
waveguide  mode.  In  general 

< h < oi/iTe  (2-19) 

In  terms  of  the  phase  velocity  w/k,  this  means  that  a wave  travels  fastest 
in  free  space  where  c = eQ.  But  a wave  propagating  along  a dielectric  guide 
moves  faster  than  a wave  traveling  through  an  infinite  dielectric.  Each  of  the 
guided  modes  propagates  with  a different  phase  velocity. 

By  renormalizing  Ey,  we  have 

Ey  = AcosBx  | x | < d 

Ey  = Acospd  e'Y(x-d)  x > d 

Ey  = AcosBd  e+v^x‘cl^  x < -d 

The  field  profile  for  the  lowest  order  mode  as  a function  of  x is  sketched 
in  Figure  2-2.  The  important  thing  to  notice  is  that  the  fields  evanesce  outside 
the  slab.  This  means  that  no  energy  is  lost  to  radiation  at  large  distances  from 
the  slab.  However,  a substantial  portion  of  the  energy  is  propagating  outside  the 
slab. 


(2-20) 
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Instead  of  choosing  fields  that  decay  with  increasing  x,  we  could  let 

f (x ) = Ae1px  + Be'1px  (2-21a) 

outside  of  the  slab,  and 

f(x)  = Csinox  ( 2-21 b ) 

inside  the  slab.  Repeating  the  analysis,  this  leads  to 

Ey  = Ae1px  + Be’1px  (2-22a) 

Ey  = Ccosax  (2-22b) 

where  A,  B,  and  C have  been  renormalized,  and 


a = (w2u0e  - h2) 


2x1/2 


P = - h2) 


2x1/2 


(2-23) 


o o 


The  propagation  constant  p can  assume  any  value.  Unlike  the  previous 
example  where  h is  a discrete  eigenvalue  determined  by  (2-18),  h can  now  have  any 
value  in  the  range 


0 < h <_  o)/m0e0  (2-24) 

In  general,  there  are  two  sets  of  solutions  for  the  fields  of  a dielectric 
waveguide.  The  first  set  includes  all  guided,  or  bound  modes,  which  are  a finite 
set  of  waves  characterized  by  discrete  values  of  the  propagation  constant  h. 

These  modes  evanesce  outside  the  guide  and  do  not  radiate.  The  constant  h is 
determined  by  a transcendental  equation  such  as  (2-18). 

The  second  set  of  solutions  is  comprised  of  the  radiation  modes,  or  modes  of 
the  continuum.^  They  are  infinite  in  number,  and  can  be  described  as  incident 
waves  that  are  reflected  at  the  dielectric  in  such  a way  that  all  boundary 
conditions  are  satisfied.  The  value  of  h is  now  continuously  variable,  and  does 
not  depend  on  a transcendental  equation. 

11.  Marcuse,  D.,  "Mode  Conversion  Caused  by  Surface  Imperfections  of  a Dielectric 
Slab  Waveguide,"  B.S.T.J.  48,  No.  10,  pp.  3177-3215,  December  1969. 
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The  guided  and  radiation  modes  are  analogous  to  the  quantum  mechanical 
solution  for  the  wavefunctions  of  an  electron  in  a potential  well.  The  guided 
modes  or  radiation  modes,  correspond  to  free  electrons,  which  can  have  any  energy. 

This  analogy  isn't  as  farfetched  as  it  seems.  For  example,  it  can  be  shown 

that  the  angular  momentum  per  unit  energy  for  a guided  mode  on  a dielectric 

cylinder  is  quantized  in  units  of  -nh',  where  n is  an  integer  and  h'  is  Planck's 

12 

constant  divided  by  2n.  However,  the  subject  of  quantized  fields  is  not 
essential  to  the  understanding  of  dielectric  waveguide  and  will  not  be  considered 
here. 

Under  appropriate  conditions,  the  guided  and  radiation  modes  of  the  dielectric 
waveguide  can  couple.  However,  for  a perfectly  uniform  and  lossless  guide  with  no 
bends  or  surface  imperfections,  the  modes  will  not  interact. 

2 . 2 THE  DIELECTRIC  CYLINDE R 

The  boundary  value  problem  of  the  cylindrical  guide  has  been  solved  by 
several  authors.  The  derivation  presented  here  follows  Stratton,  and  utilizes 
the  same  approach  that  was  used  for  the  dielectric  slab  of  section  2.1 . ^ 0 
However,  here  we  will  find  a guidance  condition  that  describes  all  of  the  modes, 
rather  than  limiting  the  argument  to  a specific  type  of  mode  (i.e.,  TE  or  TM). 

Consider  the  dielectric  cylinder  of  Figure  2-3.  The  cylinder  has  dielectric 
constant  ei  and  permeability  pi . The  surrounding  medium  is  characterized  by 
e = c2  and  u = y2.  The  radius  of  the  cylinder  is  a. 


We  can  define  a Hertz  vector  n = II*  = ipz , where 


if.  = f(r,0)e±lhz  ~ iwt 

Equations  (2-5a)  and  (2-5b)  are  now  identical  and  become 


(2-25) 


3..  Tr  If  1 + I 

3r  3r  r 30> 


+ (k2  - h2)f  = 0 


(2-26) 


Using  a separation  of  variables  technique,  let  f(r,6)  = fi(r)f2(0).  Then 


10.  Stratton,  J.  A.,  Electromagnetic  Theory,  McGraw-Hill,  New  York,  1941. 
12.  Kapany,  N.S.,  and  Burke,  J.J.,  Optical  Waveguides,  Academic  Press,  New 
York,  1972,  p.  96. 
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i (. "-)  • (. 

dr  \ dr  / \ 


(k2  - h2)r2  - p2 


)f" 


(2-27) 


+ p*f2  = 0 


(2-28) 


where  p is  a separation  constant.  Equation  (2-28)  is  Bessel's  equation,  and  has 
solutions 


Mr)  = Zp(/  k2  - h2r) 


(2-29) 


where  Zp  can  be  any  of  the  following 


Op(Xr) 

Bessel  function 

(2-30) 

Np(Xr) 

Neumann  function 

(2-31) 

HjUr) 

Hankel  function  of 
first  kind 

(2-32) 

H-(Xf) 

Hankel  function  of 

(2-33) 

p 

second  kind 

where  X = (k2  - h2)*5. 

Equation  (2-28)  is 

the  wave  equation  and  has  solutions  of  the  form 

Me)  = e"* 

(2-34) 

i 
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The  general  solution  for  t|>  can  be  expressed  as 


et"0 


(J  k2  - h2r)e 


±ihz  - iwt 


(2-35) 


where  the  constant  h is  determined  by  applying  the  boundary  conditions  at  the 
surface  of  the  cylinder  to  appropriate  solutions  of  Zn> 

From  the  relationships  (2-4)  we  have 


E1  = V x V x 

<j)Z  (2- 36a) 

H 1 = e|r  V x \pz 

a L 

(2- 36c ) 

f i 3 a 

E = -u^V 

x a-z  (2-36b) 

H 2 = V x V x iJjZ 

(2-36d) 

In  a cylindrical  coordinate  system 

, equations  (2-36)  reduce  to 

■ *«& 

pi  _ iih3i|i 

9 " r36 

E‘  = (k2  - h 2 )* 

(2-37a) 

li  i - i k 2 3^ 

Hr  " ucor90 

ui  _ i k3 

0 via)3r 

H’  = 0 

(2-37b) 

r-2  _ i V*U)34- 

r r30 

E0  = 'iwa!r 

E2  = 0 

(2-38a ) 

h2  = +Mi 

n6  “r30 

H2  = (k2  - h2  )g> 

(2-38b) 

Note  that  the  solutions  are  again  broken  up  into  TE  and  TM  solutions,  as 
they  were  for  the  slab  waveguide.  Also  note  that  3^/30  = i nd>n . In  order  for 
the  fields  to  be  single-valued  functions  of  0,  it  is  necessary  for  n to  be  an 
integer.  Since  this  is  the  case,  it  is  possible  to  write  the  general  solution 

for  the  fields  as  a summation  over  the  mode  index  n from  n = -°°  to  «.  The  mode 

amplitudes  are  an  for  the  TM  components  and  bn  for  the  TE  components.  The  general 
solutions  for  the  electric  field  are 
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Jin  . HS>  y b 

nTr  r L n n 

(2- 39a) 

"V„  -iu“Jbnf? 

( 2- 39b ) 

?an*n 

(2-39c) 

where  all  of  the  summations  are  from  n = -*>  to  The  equations  for  the  maqnetic 
field  are  similar. 

These  equations  are  the  starting  place  for  the  analysis  of  optical  fibers, 
and  are  referred  to  in  many  journal  articles. 

To  solve  for  the  fields  of  the  dielectric  cylinder,  we  must  choose  appropriate 
forms  for  in  each  of  two  regions,  r < a and  r > a.  Then  we  must  use  the  boundary 
conditions  to  match  the  solutions  at  r = a.  For  a satisfactory  solution,  the 
fields  E and  H must  decay  for  large  r,  and  be  finite  on  the  axis.  To  insure  this, 
let 

J'n  = Jn(AY-  h2r)  (2-40) 

for  r < a 


and 


il'n  = - h’r)  (2-41) 

for  r > a 

The  Bessel  function  is  finite  at  r = 0,  and  the  Hankel  function  of  the 
first  kind  decays  for  large  values  of  r.  It  is  convenient  to  keep  Stratton's 
notation  and  define 
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k = k i - u)»  u j £ i 


inside  the  cyl inder 


k - k 2 = (O*  VI 2 C ; 


outside  the  cylinder 


For  r < a,  the  z components  of  the  fields  are 


(2-42a) 


(2-42b) 


Ez  = KJn('>r>a>n 
n 


(2-43a) 


n 


(2-43b) 


and  for  r > a. 


Ez  = ^Hi(X2r)a^]Fn 


(2-44a ) 


Hz  * ZCH;(A2r)b*]Fn 


(2-44b) 


where  \\  = ki  - hJ , \*  = k|  - h2 , and 


Fp  = exp{inO  + ihz  - icot) 


(2- 45) 


The  r and  6 components  can  be  derived  in  a similar  manner. 

The  boundary  conditions  are  that  tangential  E and  H must  be  continuous  at 
r = a.  Matching  the  z components  of  E and  H yields 


Jn(uK  = Hn(v)an 


(2-46a) 
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J„(u)b’  • H‘(v)b* 


(2-46b) 


The  0 components  are  matched  in  the  same  way.  However,  the  expressions  are 
more  complicated.  We  get 


nh  . ipito  . nh  ip2u  , p 

— J„(“K  * ^(“)b„  - - H (v)>  ♦ — Hn-(v)bn 

U U V V 


(2-47a) 


ik? 


UiCJU 


. nh  . ik|  Q nh 

"’"K  = — h;-(v)^  - - H-(v)b' 


( 2 -47b ) 


u = Xia  v = X2a 


(2-48) 


The  presence  of  a prime  denotes  differentiation  with  respect  to  the  argument 
of  the  function. 

The  set  of  equations  (2-46)  and  (2-47)  are  a linear  system  of  equations  that 

are  satisfied  by  the  mode  coefficients  a1,  b’,  a®,  and  b®.  To  have  a nontrivial 

n n n n 

solution,  their  determinant  must  be  zero.  This  gives  the  guidance  condition  for 
the  dielectric  rod. 

Pi  J'(u) 
u Jn(u) 

where  ft  = ( 1 / v2 ) - (1/u2),  and  the  superscript  1 has  been  deleted  from  the  Hankel 
functions  of  the  first  kind. 


P2  H'( v) 


v Hn(v) 


ki  j:(u) 


k I 


«„(*) 


PiU  Jn(u) 


P2V  Hn(v) 


n~hzft 


(2-49) 
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1 


This  transcendental  equation  determines  the  mode  propagation  constant  h. 
Referring  to  equations  (2-46),  (2-47),  and  (2-49),  there  are  several  signifi- 
cant observations  that  can  be  made.  If  either  the  set  of  coefficients  a„  or  b 

n n 

are  arbitrarily  ;et  to  zero,  then  the  equations  have  no  solution  except  for  the 

case  n = 0.  This  means  that  for  n / 0,  only  hybrid  solutions  with  both  Ez  / 0 

and  Hz  f 0 are  allowed.  Because  of  this  hybrid  nature,  modes  are  described  by 

the  notation  EH  and  HE  , depending  on  whether  TE  or  TM  waves  dominate  the 
12  nm  nm  3 

solution.  The  subscript  n refers  to  the  order  of  the  wavefunction.  For  each 
value  of  n,  there  is  an  infinite  number  of  roots  for  the  propagation  constant  h. 
These  roots  are  designated  by  the  subscript  m,  where  m is  related  to  the  number 
of  nulls  in  the  field  pattern. 

For  n = 0,  the  guidance  condition  reduces  to  two  independent  equations.  These 


Ui  J,(u)  U 2 Hi (v) 
u JQ(u)  v Hj(v) 


(2-50) 


Ji  (u) 


H}(v) 


Mi u Jq(u)  u2v  Hq(v) 


(2-51) 


where  the  relationships  Jq(u ) = Ji(u)  and  Hg(v)  = Hi(v)  have  been  used. 

From  equations  (2-47)  with  n = 0,  it  can  be  seen  that  ag  and  bg,  the  mode 
coefficients  for  the  n = 0 mode,  are  independent  of  each  other.  In  addition,  the 
fields  have  no  0 dependence.  Equation  (2-50)  describes  the  EHgm  modes,  which  are 
transverse  electric  and  are  also  designated  as  TEgm.  Equation  (2-51)  describes 
the  HEgm  modes,  which  are  transverse  magnetic  and  are  designated  as  TMQm.  Unlike 
the  higher  modes,  which  are  hybrid  in  nature,  the  TEgm  and  TMgm  modes  can  exist 
independently  of  each  other. 


13.  Ibid,  p.  118. 
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The  simplest  n = 1 modes  are  the  EH^  and  HE^  modes,  which  for  our  purposes 
can  be  considered  identical . These  modes  have  a unique  property  that  can  be  seen 
by  investigating  the  behavior  of  the  propagation  constant  h as  the  rod  diameter 
is  decreased. 

For  a dielectric  rod  that  is  infinite  in  diameter,  all  of  the  modes  will 
propagate.  As  the  diameter  is  decreased,  the  values  of  h corresponding  to  higher 
order  modes  become  imaginary,  and  the  modes  are  cut  off.  For  a sufficiently  small 
radius,  only  the  TEg-| , TMg-j , and  HE^  modes  will  propagate.  There  is  a critical 
radius  below  which  even  the  TEg^  and  TMg^  modes  will  be  cut  off.  However,  the 
HE^  mode  has  no  cutoff,  and  in  theory  can  propagate  on  an  infinitesimally  small 
rod.  For  this  reason,  the  HE^  mode  is  called  the  dominant  mode.  It  is  also 
referred  to  sometimes  as  the  dipole  mode. 

The  critical  radius  below  which  all  modes  except  the  HE.,  mode  are  cut  off 

14 

is  given  by 

J0(x)  = 0 (2-52a) 

x = 2ira(er  - l^Ag  ( 2 - 52b ) 

where  a is  the  critical  radius 

\g  is  the  free  space  wavelength  of 
the  propagating  wave 

e is  the  relative  dielectric  constant  €?/e\ 

In  this  report,  we  will  consider  only  dielectric  rods  smaller  than  this 
critical  diameter.  We  will  be  concerned  only  with  the  propagation  characteristics 
of  the  dominant  HE^  mode.  The  field  pattern  of  the  HE^  mode  is  sketched  in 


I 


i 

i 


i 

! 
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The  experimental  analysis  will  be  based  on  a dielectric  rod  suspended  in  air. 
Thus  we  can  let 


Ml  = U2  = 

(2-53 a) 

E2  = Eg 

(2-53b) 

k2  = CL)(Pg£g)'S  = 27T/Xg 

(2-53c) 

ki  = G}(uiEi  )'2  = 2tt /X i 

(2-53d) 

where  Xg  is  the  free  space  wavelength,  and  Xi  is  the  wavelength  we  would  expect 
if  the  dielectric  medium  ei  extended  infinitely  in  all  directions. 

It  is  convenient  to  let  h,  the  propagation  constant  along  the  axis  of  the  rod, 
be  defined  as 


h = 2t t/X  (2-53e) 

where  X without  a subscript  refers  to  the  characteristic  wavelength  of  a wave 
propagating  down  the  rod. 

The  guidance  condition  for  the  HE^  mode  is  plotted  in  Figure  2-5  as  a 
function  of  X/Xg  and  rod  diameter.  It  is  plotted  for  three  values  of  dielectric 

constant  corresponding  to  fused  quartz,  polystyrene,  and  teflon.  From  the  figure, 

u 

it  can  be  seen  that  the  phase  velocity  approaches  the  limiting  value  of  c/(e  ), 
where  er  = e/Cg,  and  c is  the  speed  of  light. 

We  have  not  solved  for  the  radiation  modes  of  the  dielectric  cylinder.  These 
modes  are  derived  from  solutions  for  that  do  not  decay  for  large  r (such  as  the 
Neumann  functions).  However,  instead  of  being  plane  waves,  as  was  the  case  for  the 
dielectric  slab,  the  radiation  modes  will  be  a superposition  of  cylindrical  waves. 
The  radiation  modes  of  the  dielectric  cylinder  are  discussed  in  more  detail  by 
Marcuse. ^ ^ 

TIT  Marcuse,  D. , "Radiation  Losses  of  the  Dominant  Mode  in  Round  Dielectric 
Waveguides,"  B.S.T.J.,  pp.  1665-1693,  October  1970. 
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Fig.  2-5  Propagation  Characteristics  of 
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2.3  ATTENUATION 

The  theory  of  the  previous  sections  assumed  that  the  dielectric  material  and 
surrounding  medium  were  lossless.  This  is  certainly  not  the  case  for  a real 
structure.  An  actual  dielectric  guide  will  have  absorption,  scattering,  and 
radiation  losses.  If  the  rod  diameter  is  above  the  cutoff  diameter,  there  will 
also  be  some  energy  exchanged  between  the  different  modes,  resulting  in  mode 
conversion  loss. 

To  make  an  experimental  evaluation  of  dielectric  waveguide,  it  will  be 
necessary  t.o  measure  the  attenuation,  and  where  possible,  to  determine  its  source. 
To  model  the  loss  theoretically,  we  will  treat  the  effects  of  loss  as  a perturba- 
tion to  the  propagation  characteristics  of  the  ME,^  mode.  Specifically,  we  will 
assume  that  the  mode  propagation  constant  h is  unchanged.  All  forms  of  loss  will 
be  lumped  together  as  an  equivalent  loss  tangent  of  the  material,  and  will  lead 
to  an  effective  attenuation  factor  a. 

Consider  a plane  wave  propagating  in  an  infinite  dielectric  medium  where 


The  wave  is  determined  by 

:\’E 
?t2 

E will  be  of  the  form 

E = EqOXpOyz  - imt)  (2-b6) 

where  y2  = -w2ur.  Letting  > ■ a + i B . 


1 ,v'[ 

pc  at2 


(2-S4 ) 


(2-55) 


a2  - p2  + 2 lap  ■ 


I 


tr- 


36 


-io2pr  1 t 1,o 2 pe  1 1 


(2-57) 
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Assuming  low  loss,  a <<  0 and  e"«e',  we  get 

3 a o>/ye  1 2a  a w/ye  ' (e ' ' /e 1 ) (2-58) 

Defining  e/Eq  = er,  and  e 11 /e'  = tan  <5,  and  assuming  e'  = e,  equation  (2-58) 
becomes 

2a  = tan  5 (2-59) 

oi/ygEg  is  the  free  space  wavevector,  and 

(o/UgCg  = 2tt/Xq  (2-60) 

The  final  result  is 

2a  - (2tt/X0)/e^  tan  6 (2-61 ) 

where 

E = Egexpf-az  ± i Bz  - io)t)  (2-61) 

a is  the  voltage  attenuation  per  unit  length,  and  2a  is  the  power  attenuation 
per  unit  length. 

We  can  relate  this  to  the  round  dielectric  guide  by  equating  6 with  h. 

However,  the  attenuation  factor  a was  derived  by  assuming  that  the  dielectric 
medium  extended  infinitely  in  all  directions.  In  the  case  of  a dielectric  guide, 
part  of  the  power  is  propagating  outside  of  the  guide  and  won't  be  directly 
attenuated  by  the  dielectric.  But  since  the  field  strength  outside  the  guide  is 
determined  by  the  field  strength  inside  the  guide,  power  loss  within  the  dielectric 
will  result  in  attenuation  of  both  the  internal  and  external  fields.  One  can  take 
this  into  account  by  applying  Poynting's  Theorem.  Consider  the  section  of  wave- 
guide in  Figure  2-6.  If  all  of  the  power  is  propagating  in  the  direction  of  the 
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dA-j  = f>dfidOz 

E«H,  dA  = — E()H/;) 

E • J - r(E^Ep  + £()E()  + EzEf ) 


FIG.  2-6  SECTION  OF  DIELECTRIC  WAVEGUIDE 
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rod,  then  the  net  power  flow  into  the  volume  element  dV  is 


Pin  ' fC  * fl  ’ dfll 


(2-62a) 


The  power  leaving  the  volume  is 


Pout  - / E * H • dS. 


(2-62b) 


The  energy  dissipated  is 


/e  • 


J dV,  where  J = aE 


(2-62c) 


(o  need  not  be  an  Ohmic  conductivity;  in  the  case  of  a lossy  dielectric, 

a = ioe  tan  6 . ) 

Conservation  of  power  requires  that 


/ E x H • dA]  - fl  x H • dA2  = fl  • J 


dV 


(2-62d] 


Defining  4>  = ft  x H • dA,  and  taking  the  limit  as  Az  -+  0,  equation  2-62d 
becomes 


a v = / e • J dA 

Redefining  a to  be  the  power  attenuation  factor  (as  opposed  to  voltage 
attenuation),  we  have 


(2-62e) 


1 d$ 
♦ dz 


(2-63a) 


This  leaves 


ot  = 


fl  • J dA 
fl  x H • dA' 


2tt  w 

f J E • J pdpd0 
0 0 

\ l Vo  - V,  ^ 


(2-63b ) 


11  


[-  I 


j 
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If  the  medium  outside  the  dielectric  is  lossless,  the  current  density  J 
will  be  zero  for  p >a.  Thus  the  integral  in  the  numerator  need  only  be  evaluated 
between  0 < p < a. 

By  introducing  the  free  space  impedance  Hq  = /pq/eq,  and  relating  the  con- 
ductivity o to  the  loss  tangent  of  the  dielectric,  equation  2-63b  can  be  rewritten 


2tt 


a = opR  = er  tan  5 WUqF.qR  = ~ cr  tan  6 R 


(2-64a) 


where 


f I |f|2/n0Pdpde 


R = 


2tt 


(2-64b) 


f J pdpd0 

Jq  Jq  P 0 0 P 


R can  be  calculated  for  the  HE i -j  mode  from  the  following  equations 


e - 1 f2  + (1/p2)  - (l/p”)  4UV 

+ (U2  + V2)X  + 


16 


R = 


(1/P2)  + 0/q2) 


2V 


2V 


UX(e  + V2)  + UY ( 1 + V2)  + — (e  + U2)  (1  + U2 ) 


(2-65) 


((p2  + q2)/(c  - ]))h  = (2TTa/XQ) 

(2-66) 

u = Xq/X  = ((eq2  + p2 )/  (q2  + p2)}^ 

(2-67) 

V = C(cf  + g)/(f  + g))'s 

(2-68) 

X = f2  + (2f  + l)/p2  - 1/p4 

(2-69) 

16.  Elsasser,  W.  M. , "Attenuation  in  a Dielectric  Circular  Rod,"  J.  Appl . Rhys.  20, 
p.  1193,  December  1949. 
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Y = -g2  - (2g  - 1 )/q2  + 1/q1* 


(2-70) 


Jn(P) 


pJn(p) 


Hj'(q) 


qHj(q) 


(2-71) 


The  significance  of  equations  (2-65)  through  (2-71)  is  that  we  can  relate 
the  observed  attenuation  factor  of  a dielectric  rod  to  the  bulk  loss  tangent  of 
the  dielectric  material. 

In  practice,  R is  found  by  measuring  the  rod  diameter,  the  guide  wavelength 
A,  and  the  free  space  wavelength  AQ.  This  allows  one  to  calculate  e from  the 
guidance  condition  (2-49),  after  which  the  remaining  equations  can  be  solved. 

The  attenuation  in  a dielectric  cylinder  is  directly  proportional  to  R.  R as  a 
function  of  rod  diameter  is  plotted  in  Figure  2-7  for  polystyrene  with  a loss 
tangent  of  .001. 

2 . 4 RADIATION  LOSS 

Although  it  is  a critical  issue,  the  problem  of  radiation  from  bends  and 
surface  imperfections  is  very  complicated,  and  only  a basic  outline  of  the  problem 
will  be  presented  here.  This  disucssion  is  based  on  the  work  of  Marcatili  and 
Marcuse  at  Bell  Laboratories. ^’^>8,9  are  pr nc-f pa] -jy  concerned  with  optical 

fibers,  but  their  theories  are  general  and  apply  at  millimeter  wave  frequencies. 

In  fact,  Marcuse  and  Marcatili  have  verified  some  of  their  ideas  experimentally 
using  teflon  rods  at  55  GHz.^ 

6.  Marcatili,  E.  A.  J.,  "Bends  in  Optical  Dielectric  Guides,"  Bell  System 
Technical  Journal,  pp.  2103-2132,  September  1969. 

7.  Marcuse,  D.,  and  Derosier,  R.M.,  "Mode  Conversion  Caused  by  Diameter  Changes 
of  a Round  Dielectric  Waveguide,"  B.S.T.J.  48,  No.  10,  pp.  3217-3232, 
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Radiation  can  be  caused  either  by  surface  imperfections  or  bends  in  the 
dielectric  rod.  Mathematically,  a bend  is  just  the  limiting  case  of  a very  long 
surface  imperfection.  The  mechanism  for  radiation  is  the  coupling  of  the  discrete, 
guided  modes  to  the  continuous,  unguided  modes.  This  occurs  at  surface  imperfec- 
tions and  bends  due  to  the  impossibility  of  matching  all  of  the  boundary  conditions 
using  only  guided  modes. 

The  radiation  from  bends  can  be  understood  by  considering  the  behavior  of  an 
equiphase  wave  front  of  a guided  wave  going  around  a bend.  The  fields  extend 
outside  of  the  rod,  as  shown  in  Figure  2-8.  Much  like  skaters  playing  crack  the 
whip,  the  part  of  the  wave  furthest  outside  the  rod  must  increase  its  velocity 
the  most  to  maintain  the  phase  of  the  wave  front.  To  maintain  equiphase,  we  must 
have 


hRA0  > k(R  + xr)A0  (2-72) 

where  h is  the  propagation  constant  of  a guided  wave  and  k is  the  propagation 

constant  of  a wave  in  free  space.  R is  the  radius  of  curvature  of  the  rod,  and 

x is  the  radial  distance  from  the  center  of  the  rod. 
r 

At  the  position 

xr  = (h  - k ) R/k  (2-73) 

the  velocity  of  the  wave  front  is  equal  to  the  speed  of  light.  To  make  it  around 
the  curve,  all  of  the  energy  in  the  region  x > xr  must  propagate  faster  than  the 
speed  of  light.  This  is  impossible,  and  all  of  the  energy  in  this  region  is 
radiated.  Marcatili  and  Miller  have  proposed  an  approximate  model  to  determine  the 
amount  of  attenuation  due  to  radiation  one  can  expect  from  a bend  in  a slab  wave- 
guide.^ They  define  a characteristic  length 

zc  = a2/2X  (2-74) 

17.  Marcatili,  E.  A.  J.,  and  Miller,  S.  E. , "Improved  Relations  Describing 

Directional  Control  in  Electromagnetic  Wave  Guidance,"  B.S.T.J.,  pp.  2161- 
2188,  September  1969. 
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zc  is  the  distance  over  which  a beam  of  light  will  remain  collimated  if  it  has  a 
transverse  field  width  a and  wavelength  X.  They  relate  this  to  the  power  decay 
rate  exp(-az),  where  z is  the  longitudinal  distance  along  the  waveguide  at  the 
bend. 

The  amount  of  energy  radiated  is  E^/Ey,  where  E^  is  the  amount  of  energy 
propagating  at  x > xr,  and  Ey  is  the  total  energy.  In  a two  dimensional  system. 


E^  = F^(x)dx 

(2-75a) 

Ey  = F2(x)dz 

(2-75b) 

2 

The  quantity  F (x)  corresponds  to  the  energy  density  of  the  field  distribution. 
For  small  a,  exp(-az)  = 1 - 2az.  The  fractional  power  loss  is 

2azc  - EL/Ey  (2-7 6a ) 

or 

a = EL/(2zcEy)  (2-76b) 

It  should  be  pointed  out  that  this  a corresponds  to  power  attenuation  due  to 

radiation,  unlike  the  previous  section,  where  2a  was  the  power  attenuation  due  to 

loss  within  the  dielectric  material. 

Using  this  approach,  it  is  possible  to  estimate  the  radiation  loss  of  a 

rectangular  dielectric  waveguide  using  the  field  solutions  for  the  unbent  guide 

to  solve  for  E^  and  Ey.  A more  rigorous  approach  to  the  problem  of  radiation 

from  a rectangular  guide  has  been  used  to  verify  the  validity  of  this  approximate 
18 

method.  However,  this  method  is  not  very  useful  for  the  circular  dielectric 
guide  because  of  the  problem  of  solving  for  EL  and  Ey  in  a three  dimensional 


18.  Marcatili,  E.  A.  J.,  "Bends  in  Optical  Dielectric  Guides,"  B.S.T.J., 
pp.  2103-2132,  September  1969. 
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geometry.  It  is  also  difficult  to  determine  an  appropriate  value  for  the 
correlation  length  a.  However,  this  approach  does  offer  some  insight  into  why 
radiation  occurs. 

For  a large  diameter  rod  capable  of  propagating  more  than  one  mode,  mode 
coupling  at  a bend  will  often  be  more  important  than  radiation.  Since  we  are 
concerned  only  with  single-mode  guides,  this  isn't  a problem. 

More  important  than  bends  is  the  problem  of  radiation  from  surface  imper- 
fections. It's  possible  to  make  a straight  dielectric  guide,  but  there  will 
always  be  some  surface  roughness,  no  matter  how  smooth  the  finish.  This  problem 

has  been  studied  extensively,  and  involves  the  solution  of  coupled  mode 
19  20  21 

equations.  ’ ’ The  theory  is  straightforward.  One  solves  the  boundary  value 

problem  using  an  infinite  sum  of  guided  and  unguided  modes  to  match  boundary 
conditions  at  the  irregular  surface.  If  a single  guided  mode  is  originally 
propagating  on  the  line,  the  degree  to  which  it  couples  to  other  modes  because 
of  the  surface  imperfections  can  be  calculated  using  Fourier  techniques.  The 
mathematics  of  this  is  beyond  the  scope  of  this  report,  but  we  will  quote  two 
significant  results. 

The  surface  imperfections  on  a dielectric  guide  can  be  described  mathematically 

22 

as  a wall  distortion  function.  The  function  is  just  a mathematical  representa- 
tion of  the  surface  profile  of  the  guide.  Consider  the  slab  waveguide  of  Figure 
2-9.  If  the  height  of  the  surface  S above  the  z axi's  is  f(z),  and  the  average 
thickness  of  the  slab  is  d,  the  wall  distortion  function  is  defined  by 

W = f ( z ) - d (2-77) 


19.  Snyder,  A.  W. , "Coupling  of  Modes  on  a Tapered  Dielectric  Cylinder," 
I.E.E.E.  Trans.  Microwave  Theory  and  Techniques,  18,  No.  7,  pp.  383-392, 
July  1970. 

20.  Collin,  R.  E.,  Field  Theory  of  Guided  Waves,  McGraw-Hill,  New  York,  1960. 

21.  Marcuse,  D. , "Mode  Conversion  "Caused  by  Surface  Imperfections  of  a 
Dielectric  Slab  Waveguide,"  B.S.T.J.  48,  No.  10,  pp.  3177-3215,  December 
1969. 

22.  Marcuse,  D. , "Radiation  Losses  of  Dielectric  Waveguides  in  Terms  of  the 
Power  Spectrum  of  the  Wall  Distortion  Function,"  B.S.T.J.  48,  No.  10, 
pp.  3233-3242,  December  1969. 
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The  Fourier  transform  of  the  wall  distortion  function  is 

L 

[f (2 ) - d]e'iezdz  (2-78) 

0 

where  L is  the  length  of  the  slab,  and  0 represents  a spatial  frequency.  The 
spatial  frequencies  9 will  cause  radiation  loss  only  if 

h - k < 9 < h + k (2-79) 

where  h is  the  propagation  constant  of  the  mode,  and  k is  the  free  space  wave- 
vector  k = wi'UqEq. 

If  the  Fourier  transform  of  the  wall  distortion  function  has  no  spatial 

frequencies  in  this  range,  then  there  will  be  no  radiation  loss.  This  result 

is  not  as  simple  as  it  seems,  and  an  interested  reader  is  urged  to  consult  the 

22 

original  literature.  The  important  conclusion,  however,  is  that  surface 
imperfections  don't  necessarily  result  in  radiation,  and  that  radiation  may 
turn  out  not  to  be  a problem. 

Marcuse  has  also  calculated  the  radiation  losses  from  a round  dielectric 

23 

guide  propagating  the  HE^  mode.  He  has  included  the  effects  of  random  wall 
imperfections  in  terms  of  a correlation  function.  Let  the  power  propagating  on 
a dielectric  rod  be  P.  The  power  attenuation  constant  is  related  to  the  frac- 
tional power  loss  AP/P  by 

a = (1/L)AP/P  (2-80) 

where  L is  the  length  of  the  rod. 

It  is  possible  to  find  the  normalized  power  loss  as  a function  of  the  RMS 
deviation  in  the  radius  of  the  rod  and  a correlation  length.  If  is  the  RMS 


22.  Marcuse,  D. , "Radiation  Losses  of  Dielectric  Waveguides  in  Terms  of  the 
Power  Spectrum  of  the  Wall  Distortion  Function,"  B.S.T.J.  48,  No.  10, 
pp.  3233-3242,  December  1969. 

23.  Marcuse,  D. , "Radiation  Losses  of  the  Dominant  Mode  in  Round  Dielectric 
Waveguides,"  B.S.T.J.,  pp.  1665-1693,  October  1970. 
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deviation  in  the  radius  a of  the  rod,  the  normalized  power  loss  is 

P = (a3/A2L) (AP/P) 


(2-81) 


Marcuse  has  calculated  P as  a function  of  B/a  for  a teflon  rod  using  an 
exponential  correlation  function,  where  B is  the  correlation  length.  This  is 
reproduced  in  Figure  2-10  for  several  values  of  ka,  where  k = 2tt/Xq.  The 
correlation  length  B for  a given  guide  will  generally  be  unknown.  But  the 
function  P turns  out  to  be  bounded  from  above,  and  it  is  possible  to  compute 
the  worst  case  loss  for  a given  value  of  A or  vice  versa. 

As  an  example,  suppose  we  want  a teflon  guide  with  no  more  than  10  dB/km 
radiation  loss.  Let  the  radius  of  the  guide  be  1mm,  and  assume  the  worst  case 
value  of  P = .25,  corresponding  to  ka  = 1.5.  From  equations  (2-80)  and  (2-81), 
we  have 

(a3/A2L ) (AP/P)  = (a3a)/A2  = .25 
ci  = 10  dB/km  = 2.3  x 10"^  neper/meter 


a3  = 10  meter 
A = .003  mm 

Therefore,  the  guide  radius  must  be  uniform  to  within  .3".  Hopefully,  the 
value  of  P will  be  substantially  lower  than  worst  case. 

2.5  ISOLATION 

Since  the  fields  of  the  dielectric  waveguide  extend  outside  of  the  guide, 
they  may  interact  with  the  external  environment.  For  a practical  system,  this 
will  be  very  important,  and  may  pose  serious  problems. 

One  possible  solution  is  to  clad  the  dielectric  with  another  dielectric 
material,  as  is  done  with  optical  fibers.  The  dielectrJc  constant  of  the  cladding 
is  lower  than  that  of  the  actual  guide,  but  higher  thar,  the  dielectric  constant 
of  free  space.  The  field  solutions  inside  the  guide  are  basically  unchanged, 
but  the  external  fields  are  localized  by  the  cladding,  and  are  less  likely  to 
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FIG.  2-10  NORMALIZED  RADIATION  LOSS  DUE  TO  SURFACE  IMPERFECTIONS 

OF  A CIRCULAR  DIELECTRIC  GUIDE  AS  A FUNCTION  OF  CORRELATION 
LENGTH  FOR  c - 2.05 
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interact  with  the  external  surroundings.  Cladding  may  not  he  the  answer,  however, 
and  the  specific  solution  to  the  problem  will  depend  on  the  particular  application 
for  which  the  dielectric  guide  is  being  used. 

The  external  field  energy  has  been  computed  for  the  HEjj  mode.  The  total 
stored  energy  within  a given  radius  r is 

2n  r 

f f i f’F*  + iij.H'H*  r drdO 
0 0 1 

“ ■ -S — - ■ ' (2-8’) 

; / i + pnH*H*  r drdi’ 

0 0 

where  e = for  r < a,  and  c = for  r a. 

W has  been  computed  for  r = a,  2a,  ba,  and  10a  for  polystyrene,  fused  quartz, 
and  teflon.  These  calculations  are  included  in  the  appendix,  and  are  plotted  for 
r * a in  Figure  2-11.  Note  the  similarity  between  W and  R from  Figure  2-7.  W 
Is  the  stored  energy,  while  R relates  to  the  Poynting  power.  There  is  approxi- 
mately a factor  of  2 difference  between  the  two  quantities. 

2 . b POWER  HANDLING  CAPABILITIES 

The  power  handling  capabilities  of  a dielectric  waveguide  may  depend  on 
several  different  factors.  Excessive  power  absorption  may  cause  the  material  to 
increase  in  temperature  and  melt.  Or  the  dielectric  strength  of  the  material 
can  be  exceeded,  causing  a localized  breakdown.  In  an  optical  fiber,  the  wave- 
length of  the  light  is  so  small  that  scattering  is  important.  The  power  trans- 
mission is  then  limited  by  the  threshold  power  level  at  which  stimulated  Raman 
and  Brlllouln  scattering  occur. 

Ihe  dielectric  strength  of  a material  is  typically  specified  at  low  frequencies 
or  D.C.  At  microwave  frequencies,  the  dielectric  strength  may  he  substantially 
different,  causing  the  material  to  break  down  at  much  lower  field  intensities. 
Whether  this  is  indeed  the  case  for  millimeter  waves  will  be  difficult  to 
determine  until  high  power  sources  become  available. 
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The  problem  of  dielectric  heating  can  be  solved  approximately  to  yield  a 
rough  estimate  of  the  power  that  can  be  dissipated  by  a cylindrical  dielectric 
guide.  The  approach  used  here  will  be  to  find  the  power  level  that  will  cause 
the  temperature  at  the  center  of  the  rod  to  reach  its  melting  point.  The 
electromagnetic  energy  absorbed  in  the  rod  will  be  treated  as  a uniform  source 
of  heat  throughout  a cross-section  of  the  rod.  The  heat  diffuses  to  the  outside 
of  the  rod,  where  it  is  dissipated  by  radiation  and  convection. 

Heat  conduction  in  the  steady-state  is  described  by 

VZT  = -g/ic  (2-83) 

where  T is  the  temperature,  q is  the  heat  generated  per  unit  volume,  and  k is  the 
thermal  conductivity  of  the  dielectric.  In  a differential  length  Az,  the  heat 
generated  is 


AP  = qAV  = aPAz  (2-84) 

where  a is  the  power  attenuation  per  unit  length,  AV  is  the  differential  volume 
element  na2Az,  a is  the  radius  of  the  dielectric  rod,  and  P is  the  amount  of 
power  propagating  down  the  guide.  Equation  (2-83)  becomes 

V2T  = -cxP/(ira2k)  (2-85) 


We  would  like  to  solve  this  equation  to  find  P as  a function  of  T^,  the 
temperature  at  the  center  of  the  guide.  The  maximum  value  of  P will  be  determined 
by  setting  TQ  = T , the  melting  temperature  of  the  dielectric.  This  is  one  of 
the  two  boundary  conditions  necessary  to  solve  the  problem. 

For  the  steady-state,  we  can  assume  that  all  of  the  generated  heat  is 
diffused  to  the  surface  of  the  rod.  So  at  the  surface  of  the  guide,  we  have  the 
second  boundary  condition. 


qAV 


3T 

-K  — - 

3r 


a 


(2-88) 


where  AS  is  a unit  element  of  surface  area. 
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By  symmetry,  the  0 directed  component  of  the  temperature  gradient  is  zero. 
Thermal  diffusion  along  the  axis  of  the  guide  will  be  unimportant,  so  the  z 
component  can  also  be  ignored.  Equation  (2-83)  reduces  to 


32T 

13T 

4.  

aP 

■f  - 

3r2 

r3r 

iTa2K 

(2-87) 


where  T = T«  at  r = 0,  and  = -aP/ (2-rraic)  at  r = a. 
u dr 

The  solution  for  TQ  is 

T0  = Tw  + «P/(4ttk)  (2-88) 

T is  the  surface  temperature  of  the  guide,  and  is  determined  by  the  ambient 
w 

temperature  and  the  effecti veness  of  convection  and  thermal  radiation.  The  problem 

of  finding  T is  discussed  in  detail  by  Holman,  and  only  the  results  will  be 
w 24 

presented  here.  If  T is  the  ambient  temperature  outside  the  rod,  then 


aP  a2  a 

na2  4k  2hj. 


+ T 


oo 


(2-89) 


where  h^  is  an  equivalent  heat  conduction  coefficient  that  combines  the  effects 
of  radiation  and  convection. 

We  can  define  h^  by 

ht  ' V + hc  '2-90> 

h is  determined  by  the  Stefan-Bol tzman  law  of  black-body  radiation,  and  hc 
is  derived  using  fluid  mechanics. 


24.  Holman,  J.  P.,  Heat  Transfer,  Fourth  Edition,  McGraw-Hill,  New  York,  1976. 
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+ T’)(TW  + TJ  ( 2-91  a ) 

hc  - 1.32{(Tw  - Tj/(2a)}1/4  (2-91b) 

f,  is  the  emissivity  of  the  dielectric,  and  a = 5.67  x 10~8  Watt/m2  °K4. 

Unfortunately,  h^,  hc,  and  k are  all  functions  of  temperature.  However,  the 
dielectric  rods  are  only  a few  millimeters  in  diameter  at  most,  and 


Jn  » T 
0 w 


(2-92) 


As  an  example,  consider  a 2 mm  teflon  rod  with  a loss  tangent  of  5 x 10  . 

From  the  theory  of  section  2.3,  a is  found  to  be  .0048  neper/cm.  The  melting 

point  of  teflon  is  321°C,  and  the  thermal  conductivity  is  [4.86  x 10~8T  + 2.53 
-3  25 

x 10  ]Watt/cm  °K.  Depending  on  what  value  one  assumes  for  £,  P will  be  in  the 
neighborhood  of  100  to  200  Watts.  Reducing  the  loss  tangent  of  the  material  will 
substantially  increase  this  amount,  as  will  decreasing  the  rod  diameter. 


25.  Brandrup  and  Immergut.,  editors,  PoJymer_  Handbook , Second  Edition, 
Wiley  Sons,  New  York,  1975. 
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Chapter  III 

EXPERIMENTAL  OBSERVATIONS 


3.1  INTRODUCTION 

Before  we  can  legitimately  use  the  theory  of  the  HE^  mode  to  describe  the 

behavior  of  an  actual  guide,  it  will  be  necessary  to  verify  that  the  HE^  mode 

actually  propagates  when  a wave  is  launched  on  a dielectric  rod.  From  Chapter  2, 

we  know  that  for  a given  d/Xg  and  c,  the  HE^  mode  will  be  characterized  by  a 

certain  X/Xg.  If  we  assume  that  the  HE^  mode  is  propagating,  then  measurement 

of  d/Xg  and  X/Xg  will  yield  a value  for  c.  If  this  c agrees  with  the  known  e 

for  the  material,  we  can  be  sure  that  the  HE^  mode  is  propagating.  We  can  look 

for  the  presence  of  other  modes  by  probing  the  standing  wave  pattern  on  the  rod. 

Since  each  mode  has  its  own  characteristic  wavelength,  the  presence  of  more  than 

one  mode  will  be  readily  apparent.  This  method  has  been  used  by  Cullen  and 

26 

Gillespie  to  study  the  properties  of  a glass  rod  at  33  GHz. 

Unfortunately,  values  for  e at  70  GHz  are  not  available  for  most  materials. 
However,  dielectric  constant  varies  slowly  with  frequency,  and  we  can  approximate 
e with  a value  measured  at  a lower  frequency  (i.e.,  25  GHz).  The  difference  in 
X/Xg  between  each  of  the  lowest  order  modes  is  substantial  enough  that  only 
approximate  agreement  between  the  observed  and  known  values  for  c is  necessary 
to  confirm  the  presence  of  the  HE^  mode.  Once  this  is  verified,  we  can  assume 
that  the  observed  value  of  c is  the  actual  dielectric  constant  of  the  material 
at  70  GHz.  Thus  measurement  of  d/Xg  and  X/Xg  will  yield  a previously  unknown 
value  for  e at  70  GHz,  as  well  as  verify  the  presence  of  the  dipole  mode. 

The  attenuation  characteristics  of  the  dielectric  rod  are  of  critical 
importance  in  determining  the  practicality  of  using  the  rod  as  a waveguide.  In 
Chapter  2,  it  was  seen  that  knowledge  of  the  loss  tangent  of  the  dielectric  will 

26.  Cullen,  A.  L.,  and  Gillespie,  E.  F.  F.,  "A  New  Method  for  Dielectric 
Measurements  at  Millimeter  Wavelengths,"  Presented  at  the  Symposium  on 
Millimeter  Waves,  Polytechnic  Institute  of  Brooklyn,  April  1959. 
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yield  an  attenuation  factor  for  the  guide.  However,  this  does  not  include  the 
attenuation  due  to  radiation,  which  may  be  significant.  In  most  cases,  data 
for  tan  5 at  70  GHz  isn't  available.  Thus  it  is  desirable  to  measure  the  attenua- 
tion directly.  We  can  relate  this  measured  attenuation  to  an  equivalent  tan  6' 
that  will  include  the  effects  of  radiation.  In  those  cases  where  the  intrinsic 
loss  tangent  of  the  material  is  known,  comparison  of  tan  5 with  tan  S'  will  allow 
one  to  determine  how  much  energy  is  lost  to  radiation. 

The  method  used  to  determine  the  guide  wavelength  and  attenuation  factor  is 
relatively  simple.  A klystron  is  used  to  excite  the  HE^  mode  on  a long  section 
of  dielectric  rod.  A probe  that  can  move  up  and  down  the  length  of  the  rod  on  a 
motorized  cart  is  attached.  The  movement  of  the  probe  is  synchronized  with  a 
chart  recorder  so  that  one  can  plot  the  average  power  propagating  on  the  rod  as 
a function  of  position.  The  power  decays  exponentially  as  the  probe  moves  down 
the  line,  thus  yielding  the  attenuation  factor.  If  a metal  reflector  is  attached 
to  the  end  of  the  dielectric  rod,  a standing  wave  is  set  up  and  the  guide  wave- 
length can  also  be  determined. 

Figure  (3-1)  shows  the  apparatus  used.  The  dielectric  rod  is  butted  against 
the  metal  waveguide  and  is  held  in  place  by  four  set  screws.  The  other  end  of  the 
rod  is  pulled  taught  and  held  with  a clamp.  For  the  first  experiments,  the  probe 
consisted  of  a short  piece  of  dielectric  rod  butted  into  a metal  waveguide  con- 
taining a diode  detector.  This  dielectric  probe  is  placed  in  proximity  to  the  rod 
being  measured.  In  later  measurements,  the  probe  was  a piece  of  metal  waveguide 
with  a hold  drilled  through  the  top.  The  dielectric  waveguide  passes  through  the 
hole,  thus  coupling  the  longitudinal  E field  of  the  HE^  mode  to  the  TE-jq  mode 
of  the  rectangular  guide.  Other  probe  configurations  were  tried,  but  this  was 
the  most  successful . 

It  is  essential  that  the  dielectric  rod  being  measured  be  as  long  as  possible. 
If  the  dielectric  rod  is  too  short,  the  attenuation  may  be  too  little  to  measure 
accurately.  In  addition,  the  VSWR  caused  by  a reflection  at  the  end  of  the  rod 
will  increase  significantly  and  make  the  attenuation  seem  lower  than  it  really 
is.  (The  probe  will  detect  the  sum  of  the  incident  ami  reflected  waves;  this  sum 
does  not  decay  exponentially,  as  examination  of  basic  transmission  line  theory 
will  show.  If  the  rod  is  sufficiently  long,  the  initial  wave  will  decay  enough 
so  that  a reflection  at  the  end  of  the  rod  will  be  insignificant). 
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EXPERIMENTAL  RESULTS 


The  very  first  experiment  was  to  determine  that  the  dielectric  rod  was 
indeed  acting  as  a waveguide.  This  is  the  case,  as  can  be  seen  from  Figure  3-2. 

The  top  curve  was  plotted  by  pointing  a transmitting  and  receiving  horn  at  each 
other  and  moving  the  two  horns  apart.  The  gain  of  the  horns  is  high  enough  that 
a standing  wave  pattern  is  set  up.  As  the  horns  move  apart,  the  amplitude  of  the 
received  signal  decreases  greatly.  The  length  of  the  top  curve  in  the  figure 
corresponds  to  about  15  cm  movement  of  the  receiving  horn.  The  free  space  wave- 
length can  be  computed  by  counting  the  number  of  oscillations. 

The  second  curve  was  made  by  placing  a polystyrene  rod  between  the  trans- 
mitting horn  and  a stationary  reflector.  A dielectric  probe  was  placed  in  proxi- 
mity to  the  rod  and  moved  along  its  length.  Note  that  the  amplitude  is  no  longer 
a strong  function  of  position.  In  addition,  the  characteristic  wavelength  has 
decreased,  corresponding  to  the  slower  phase  velocity  of  the  wave  traveling  through 
a dielectric  medium.  The  changes  in  amplitude  along  the  path  are  due  to  variation 
in  the  separation  between  the  probe  and  rod.  Since  we  are  looking  at  a standing 
wave  pattern,  adjacent  minima  are  separated  by  A/2,  not  A. 

The  next  experiment  was  to  measure  the  dielectric  constants  of  fused  quartz, 
polystyrene,  and  teflon.  The  results  are  shown  in  Figure  (3-3). 


MATERIAL 

Fused  Quartz 

Polystyrene 

Teflon 

Rod  diam. 

1.09  mm 

2.38  mm 

2.01  mm 

dc 

1 . 95  mm 

2. 62  mm 

3.24  mm 

d/A(j 

.259 

.563 

.477 

A/\0 

.980 

.811 

.947 

e 

3.75 

2.53 

2.00 

e(25  GHz)28 

3.78 

2.56 

2.08 

Figure  3-3  Measurement 

of  Dielectric  Constant  at  71 

.0  GHz 

28.  VonHIppel , A.,  Dielectric  Materials  and  Applications,  M.I.T.  Press  and  John 
Wiley  & Sons,  New  York,  1954. 
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The  guide  wavelength  was  measured  by  probing  a section  of  dielectric  rod.  XQ  was 
determined  by  measuring  the  frequency  of  the  klystron  with  a cavity  wavemeter.  The 
rod  diameter  was  measured  with  a micrometer,  e is  computed  using  a Fortran  program 
to  find  a best  fit  to  the  guidance  condition  of  the  HE^  mode.  The  quantity  dc 
corresponds  to  the  critical  diameter  above  which  higher  order  modes  begin  to 
propagate. 

Since  the  measured  and  known  values  for  e are  in  close  agreement  for  all 
three  materials,  we  can  be  sure  that  the  observed  mode  is  actually  the  HE^  mode. 

To  see  if  anything  unusual  happens  when  the  rod  diameter  exceeds  dc,  the 
standing  wave  patterns  were  probed  for  four  rods  with  d > dc-  Interference 
patterns  were  observed  in  all  four  cases.  The  results  are  shown  in  Figure  3-4. 

There  are  several  possible  explanations  for  what  is  going  on.  One  is  that 

the  TEqi  and  TM^  modes  are  also  propagating.  This  may  be  the  case  for  the 

3.24  mm  polystyrene.  But  for  the  teflon  and  2 mm  fused  quartz  rods,  the  beat 

length  is  too  long  to  be  mode  interference.  In  these  cases,  it  may  be  that  the 

rods  are  elliptical  in  cross-section,  instead  of  being  circular.  When  this 

happens,  the  HE^  mode  is  split  up  into  two  modes,  each  polarized  along  a different 

axis  of  the  ellipse.  The  two  polarizations  will  each  propagate  with  a slightly 

different  phase  velocity,  accounting  for  the  long  beat  length.  This  effect  is 

26 

discussed  in  more  detail  by  Cullen  and  Gillespie.  The  3 mm  fused  quartz  may  be 
experiencing  the  same  phenomena,  but  it  isn't  clear  from  the  short  length  of  the 
plot. 

The  final  experiment  was  to  measure  the  attenuation  characteristics  of 
various  rods.  This  was  done  by  probing  the  rods  for  a distance  of  one  to  1 1/2 
meters  along  their  length.  Figure  3-5  shows  the  average  power  as  a function  of 
length  for  1.5  meters  along  a polystyrene  rod.  This  particular  rod  is  made  of 
Rexolite  1422,  which  is  polystyrene  that  has  become  crossl inked  from  exposure  to 
radiation. 

Figure  3-6  shows  a tabulation  of  the  results  for  an  assortment  of  poly- 
styrene, Rexolite,  and  teflon  rods.  To  check  the  measurements  against  theory, 

26.  Cullen,  A.  L.,  and  Gillespie,  E.  F.  F.,  "A  New  Method  for  Dielectric 
Measurements  at  Millimeter  Wavelengths,"  Presented  at  the  Symposium  on 
Millimeter  Waves,"  Polytechnic  Institute  of  Brooklyn,  April  1959. 
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FIG  3 7 ATTENUATION  AS  A FUNCTION  OF  DIAMETER  FOR  . - 2.56 
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the  loss  tangent  is  computed  for  the  various  polystyrene  and  Rexolite  rods.  This 

was  not  very  difficult,  as  the  numerical  calculations  for  attenuation  as  a function 

of  rod  diameter  have  already  been  done  by  El sasser1 ^ for  e = 2.56  (see  section 

2.3).  The  results  are  plotted  in  Figure  3-7.  The  solid  lines  represent  the 

theoretical  calculations  for  various  values  of  tan  6.  The  experimental  data  also 

27 

includes  measurements  made  at  24  GHz  by  Chandler. 

Some  of  the  larger  rods  exceed  the  critical  diameter,  and  more  than  one  mode 
may  be  propagating.  Since  the  attenuation  of  each  mode  is  slightly  different,  the 
experimental  data  for  these  rods  may  not  represent  the  true  attenuation  of  the  HE^ 
mode.  However,  the  results  for  the  larger  rods  are  consistent  enough  to  warrant 
their  inclusion. 

The  computer  values  of  tan  & appear  to  be  self  consistent,  thus  supporting 

the  validity  of  El sasser's  theory.  However,  there  is  no  independent  data  for  tan 

6 of  polystyrene  at  70  GHz  with  which  to  compare  these  results.  There  is  a 

possibility  that  some  of  the  attenuation  is  due  to  radiation,  and  that  the  measured 

tan  5‘  is  higher  than  the  intrinsic  tan  <S  of  the  material.  Measurements  of  tan  .s 

28 

have  been  made  for  polystyrene  and  Rexolite  at  25  GHz,  and  for  polystyrene  at 
29 

400  GHz.  Although  the  data  reported  in  Von  Hippel  varies  substantial ly  from 
sample  to  sample,  it  is  clear  that  tan  <5  rises  dramatically  between  24  and  400 
GHz.  At  24  GHz,  tan  iS  is  roughly  5 x 10"^;  at  400  GHz,  tan  6 is  reported  as 

-3  _3 

6 x 10  . Our  experimental  data  indicates  that  tan  & is  approximately  10  at 

72.70  GHz.  The  indications  are  that  radiation  loss  is  not  significant. 

The  overall  conclusion  to  be  drawn  is  that  dielectric  waveguide  can  be 
designed  to  have  lower  loss  than  metal  waveguide  at  these  frequencies,  simply  by 
decreasing  the  diameter.  The  polystyrene  has  substantially  more  loss  than  the 


16.  Elsasser,  W.  M. , "Attenuation  in  a Dielectric  Circular  Rod,"  0.  Appl . 
20,  p.  1193,  December  1949. 

27.  Chandler,  C.  H. , "An  Investigation  of  Dielectric  Rod  as  Wave  Guide," 
Journal  of  Applied  Physics,  20,  p.  1188,  December  1949. 

28.  VonHippel,  A.,  Dielectric Materials  and  Applications,  M.I.T.  Press 
and  John  Wiley  & Sons,  New  York,  1954". 

29.  Breeden,  K.  H.,  and  Sheppard,  A.  P.,  "A  Note  on  the  Millimeter  and 
Submillimeter  Wave  Dielectric  Constant  and  Loss  Tangent  Value  of 
Some  Common  Materials,"  Radio  Science,  Vol . 3,  No.  2.,  February  1968. 
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teflon,  primarily  for  two  reasons.  The  teflon  has  a slightly  lower  dielectric 

constant,  resulting  in  less  energy  propagating  within  the  rod  for  a given  diameter. 

In  addition,  teflon  is  a nonpolar  polymer.  This  makes  it  relatively  inert  in  the 

presence  of  an  external  field.  Polystyrene,  on  the  other  hand,  is  extremely  polar. 

This  results  in  a strong  interaction  between  the  individual  molecules  and  the 

30 

applied  electric  field,  which  in  turn  causes  loss. 


30.  P.  Debye,  Polar  Molecules , Dover  Publications,  New  York,  NY,  1945. 
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Chapter  IV 

SUMMARY  AND  CONCLUSIONS 

In  this  report,  we  have  approached  the  subject  of  the  dielectric  waveguide 
from  both  a theoretical  and  experimental  point  of  view.  Even  though  the  discussion 
has  been  limited  to  the  dominant  mode  of  the  cylindrical  guide,  the  theory  is 
complicated  and  the  mathematics  difficult.  Although  they  seem  straightforward , 
the  experimental  technigues  involve  a great  deal  of  numerical  analysis,  and 
require  a thorough  understanding  of  basic  theory. 

Our  discussion  provides  a qualitative  understanding  of  the  basic  behavior 
of  a dielectric  waveguide.  The  experiments  have  confirmed  that  the  basic  theory 
is  sound,  and  that  the  use  of  dielectric  rods  for  making  millimeter  wave  trans- 
mission lines  is  not  unreasonable. 

To  summarize  briefly,  the  behavior  of  a dielectric  waveguide  is  determined 
by  the  electrical  and  magnetic  properties  of  the  material,  and  by  the  physical 
structure  of  the  waveguide.  For  a homogeneous,  non-magnetic  material,  the  major 
difficulty  lies  in  the  solution  of  Maxwell's  equations  subject  to  the  appropriate 
boundary  conditions.  The  most  important  property  of  any  dielectric  waveguide  is 
that  the  field  solutions  can  be  separated  into  a discrete  set  of  guided  modes, 
and  a continuous  set  of  radiation  modes.  Each  guided  mode  is  characterized  by 
an  eigenvalue  h and  a phase  velocity  co/k.  The  guided  modes  decay  at  large 
distances  from  the  guide  and  do  not  radiate  energy.  However,  a substantial 
portion  of  the  energy  of  the  guided  mode  propagates  outside  of  the  guide. 

The  attenuation  factor  of  the  guide  depends  on  the  loss  tangent  of  the 
dielectric  material,  and  on  the  extent  to  which  the  energy  of  the  mode  is 
contained  within  the  dielectric.  By  reducing  the  guide  dimensions,  the  relative 
amount  of  energy  propagating  outside  of  the  guide  is  increased,  and  the  attenuation 
factor  is  decreased.  The  limiting  factor  in  reducing  the  attenuation  is  the 
critical  guide  cross  section  below  which  a given  mode  is  cutoff  and  doesn't 
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propagate.  The  dominant  mode  of  the  circular  guide  has  no  cutoff,  and  in  principle, 
the  attenuation  of  this  mode  can  be  made  arbitrarily  small. 

Radiation  loss  is  caused  by  the  coupling  between  guided  and  radiation  modes. 
This  occurs  at  bends  in  the  guide  or  at  surface  imperfections . A perfectly  straight 
guide  with  no  surface  imperfections  will  not  radiate.  If  the  spatial  frequencies 
of  the  Fourier  transform  of  the  surface  profile  of  the  guide  are  not  in  a g: ven 
frequency  range,  then  radiation  due  to  surface  imperfections  will  be  minima';. 
Radiation  at  bends  is  due  to  the  impossibility  of  a wave  front  propagating  around 
a bend  at  a speed  greater  than  the  speed  of  light.  For  large  diameter  guides, 
mode  coupling  will  often  be  more  important  than  radiation. 

The  amount  of  power  absorbed  by  a dielectric  guide  depends  on  the  loss 
tangent,  the  homogeneity  of  the  material,  and  the  field  distributions  inside  the 
guide.  The  ability  of  the  guide  to  dissipate  energy  will  depend  on  the  size  and 
shape  of  the  guide,  the  thermal  conductivity  of  the  dielectric,  and  the  effective- 
ness of  convection  and  thermal  radiation  at  the  surface  of  the  guide.  For  the 
steady  state,  one  can  approximate  the  maximum  amount  of  energy  that  can  be 
dissipated  by  a guide  without  its  melting.  However,  inhomogeneities  in  the 
dielectric  will  probably  cause  localized  breakdown  of  the  rod  at  much  lower  power 
levels,  so  that  this  analysis  is  only  approximate. 

In  this  study,  we  have  experimentally  determined  the  dielectric  constants  of 
samples  of  fused  quartz,  polystyrene,  and  teflon  at  a frequency  of  71.0  GHz.  The 
results  are  consistent  with  values  of  ep  measured  at  other  frequencies.  The 
attenuation  factor  was  measured  for  various  dielectric  rods  at  71.7  GHz.  The 
attenuation  factors  for  polystyrene  were  related  to  an  equivalent  loss  tangent 
for  the  material  that  includes  the  effects  of  both  radiation  and  absorption  loss. 
Comparison  of  this  result  with  other  data  indicates  that  radiation  loss  probably 
is  not  significant. 

There  is  a major  problem  that  has  only  been  touched  on  in  this  report.  The 
manufacturing  techniques  have  a very  large  effect  on  the  electrical  properties  of 
a material,  and  there  is  often  a large  disparity  between  the  electrical  properties 
of  various  samples  of  the  same  material.  Unless  the  chemical  and  mechanical 
properties  of  a dielectric  sample  are  also  specified,  it  is  not  accurate  to  say 
that  the  measured  electrical  properties  represent  the  true  properties  of  anythinq 
except  the  particular  sample  used  in  the  experiment.  For  example,  the  degree  of 
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crystallinity  of  a teflon  rod  may  be  different  than  that  of  a teflon  slab,  because 
of  the  different  manufacturing  processes  involved.  Consequently,  the  complex 
dielectric  constant  of  a teflon  slab  will  be  slightly  different  than  that  of  a 
teflon  rod.  The  amount  of  impurities  will  have  a major  effect  on  the  electrical 
properties.  But  the  amount  of  impurities  also  varies  strongly  between  samples. 

Unless  a material  is  specifically  designed  for  electrical  appl ications , it 
is  unlikely  that  it  will  be  very  pure.  For  instance,  plasticizers  are  routinely 
added  to  many  polymers  (i.e.  - polystyrene)  to  improve  their  mechanical  properties. 
The  addition  of  a plasticizer  will  change  the  dielectric  constant  and  raise  the 
loss  tangent. 

The  materials  used  in  the  experiments  were  not  designed  with  millimeter  wave 
applications  in  mind,  and  there  is  no  doubt  that  their  properties  at  millimeter 
wave  frequencies  can  be  improved.  In  addition  to  decreasing  the  amount  of 
impurities  in  the  materials,  techniques  such  as  cross-linking  the  polymers  or 
casting  the  materials  in  a vacuum  instead  of  air  will  substantially  improve  the 
dielectric  and  attenuation  characteristics  of  the  materials. 

The  effect  of  temperature  on  dielectric  properties  is  also  a major  problem 
that  will  have  to  be  solved. 

The  point  is  that  the  development  and  understanding  of  materials  is  a critical 
issue  in  the  design  of  practical  dielectric  devices.  This  is  true  for  any  structure 
be  it  dielectric  waveguide,  stripline,  microstrip,  or  whatever. 

In  this  report,  we  have  shown  that  a dielectric  rod  waveguide  can  guide  a 
wave  with  an  amount  of  attenuation  that  is  lower  than  that  of  a rectangular  metal 
guide.  The  dielectric  waveguide  has  poor  isolation,  and  is  difficult  to  interface 
with  other  structures.  But  it  is  extremely  inexpensive  compared  to  other 
structures,  and  is  simple  from  a physical  point  of  view. 

Further  work  needs  to  be  done  in  several  areas.  The  theory  needs  to  be 
extended  to  hiqhp>"  order  modes  and  to  topologies  other  than  the  circular  rod. 

The  measurement  techniques  must  be  improved  until  they  are  simple,  accurate,  and 
versatile  enough  for  general  use  in  analyzing  materials.  The  study  of  materials 
is  of  critical  importance,  as  the  practicality  of  any  structure  will  be  determined 
by  the  material's  properties.  And  finally,  the  design  of  components  such  as 
antennas,  phase  shifters,  delay  lines,  and  couplers,  is  a wide  open  area. 
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Dielectric  waveguide  may  not  be  the  ultimate  waveguidinq  structure,  but  there 
is  no  doubt  that  it  will  be  useful  for  some  applications,  particularly  antennas 
and  long  distance  transmission  lines.  Ten  years  ago,  few  people  would  have 
believed  that  optical  fibers  would  be  practical  for  telephone  communications. 

But  the  widespread  use  of  such  fibers  is  just  around  the  corner.  Ten  years  from 
now,  the  use  of  dielectric  waveguide  at  millimeter  wave  frequencies  may  be  just  as 
practical . 
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Appendix  A 

NUMERICAL  COMPUTATIONS 

Solutions  to  the  guidance  condition  of  the  HE-ji  mode  have  been  calculated 
for  values  of  e corresponding  to  teflon  (2.08),  polystyrene  (2.56),  and  fused 
quartz  (3.78).  The  data  is  tabulated  on  the  following  pages. 

The  diameter  factor  refers  to  the  distance  from  the  rod  within  which  the 
fractional  internal  energy  was  computed.  The  calculation  of  the  fractional 
internal  energy  is  discussed  in  chapter  2. 

Z1  and  Z2  are  calculation  parameters  that  are  related  to  the  Hankel  and 
Bessel  functions  that  describe  the  electric  and  magnetic  fields. 
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